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ABSTRACT 

We compute the 7-ray and neutrino diffuse emission of the Galaxy on the basis of a recently proposed 
phenomenological model characterized by radially dependent cosmic-ray (CR) transport properties. 
We show how this model, designed to reproduce both Fermi -LAT 7-ray data and local CR observables, 
naturally reproduces the anomalous TeV diffuse emission observed by Milagro in the inner Galactic 
plane. Above 100 TeV our picture predicts a neutrino flux that is about five (two) times larger 
than the neutrino flux computed with conventional models in the Galactic Center region (full-sky). 
Explaining in that way up to ~ 25% of the flux measured by IceCube, we reproduce the full-sky 
IceCube spectrum adding an extra-Galactic component derived from the muonic neutrino flux in the 
northern hemisphere. We also present precise predictions for the Galactic plane region where the flux 
is dominated by the Galactic emission. 


1. INTRODUCTION. 


In recent years the IceCube collaboration has opened 
the era of neutrino astronomy and announced the 
detection of 37 extraterrest ria l neu trinos above ~ 
30 TeV (jAartsen et all 2013bllaL 1201411. Mo re recently, 
a preliminary analysis ( Aartsen et al.ll2015df> . based on 
four years of data, rose the total number of high-energy 
starting events (HESE) to 54. 

The astrophysical spectrum inferred by IceCube on the 
basis on the three-year data set was fitted by a power law 
with index T = —2.3 ± 0.3 above 60 TeV (jAartsen et all 
l2014fl , while the four- year data favor a steep er spectrum: 
T = —2.58 ± 0.25 (( Aartsen et all l2015dlH . Although 


a statistically significant departure from isotropy can- 
not be claimed yet, a recent analysis (|Ahlers et al.ll201~5 l 


showed that the angular distribution of HESE events al¬ 
lows up to 50% of the full-sky astrophysical flux to have 
a Galactic origin. Moreover, a hint of a harder spectrum 
in the no rthern Hemisphere ma y be suggested by a recent 
analysis (1 Aartsen et all 12015bf) . 
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The Galaxy, indeed, is a guaranteed source of neutrinos 
up to a fraction of PeV energies at least. 

A sizable flux may either come from freshly accelerated 
cosmic rays (CRs) undergoing hadronic scattering with 
gas clumps, or from the hadronic interactions between 
the Galactic CR sea and the diffuse gas. 

The former scenario, however, cannot explain the 
steepness of the neutrino spectrum measured by IceCube 
and is in tension with Fermi -LAT up per lim it on the cor¬ 
responding 7-ray emission dTchernin et alJ I2013jh 

In the latter, instead, if the local CR spectrum 
is assumed to be representative of the entire Galac¬ 
tic population, the computed spectrum should be si g¬ 
nifica n tly lower than the mea s ured spectrum dSteckerl 
Il979t iBerezinskv et al. 119931 lEvoli et all 120071) (see 
also (lAhlers et al.l 12015 1: The authors show that only 
~ 8% of the HESE can be accounted in that way, under 
the conventional assumption that the same CR transport 
properties hold throughout the whole Galaxy). 

However, it is conceivable that CR diffusion - due to 
a stronger star forming activity and peculiar magnetic 
field strength/geometry - behaves differently in the inner 
Galactic region. Several anomalies observed in the 7-ray 
diffuse emission support this possibility. 

We start by noting that conventional models cannot 
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Figure 1. Diffuse emission 7 -ray spectrum from the inner Galac¬ 
tic plane (| 6 | < 2°, 30° < l < 65°) computed for the refer¬ 
ence models considered in this Letter is compared with Fermi- 
LAT and Milagro data. The Milagro differential flux reported 
h ere is 17 % lower with respect to the flux reported in 2008 
( flAbdo et a l. 2008)) due to the assumpt ion of a spectral index 
of 2.4 instead of 2 .7 GHu ciitemcycr 20lJ)). The exp ected sensi- 
tivities of HAWC dAbevsekara et al. 2013) and CTA d Act is et al.l 
120111 ) are reported. The spectral components are shown for the 
KRA 7 model only. The Fermz-LAT data points refer to 5 years of 
data, within the event class ULTRACLEAN, according to Fermi 
tools v9r32p5. 

explain the large 7 -ray flux measured by the Milagro ob¬ 
servator y from the inner GP region at 15 TeV m edian 
energy (IProdanovic et al.l 120071 ; lAbdo et al.l l2008i) . In 
Fig. Q] we show how a representative conventional model, 
with si milar spectral properties as the Fermi benchmark 
model (lAckermann et all 120121) . clearly fails to repro¬ 
duce that measurement. This problem is common to 
all the models of this kind and still holds assuming - as 
done in (jAhlers et al.ll2015h - that the spectral hardening 
found by PAMELA in th e CR proton and he lium spec¬ 
tra above ~ 230 GeV/n (lAdriani et al.l [20l3h is present 
throughout the whole Galaxy. Therefore, the Milagro 
excess is still an open issue, and indeed its possible rele¬ 
vance for high-energy neutrino physics h a s often been 
point ed out (see e.g. (|Gabici et al.l 120081 iTavlor et al.l 

MD). 

An even more serious anomaly was found at 
lower energies in the Ferrm -LAT diffuse 7 -ray spec¬ 
trum (lAckermann et ahll 2012 l ): the conventional models 
systematically underestimate the measured flux in the in¬ 
ner GP region above a few GeV. A new ph enome nological 
scenario was proposed in (jGaggero et al.l 1201411 in order 
to account for the latter results: the idea is to consider 
a radial dependence for both the rigidity scaling index S 
of the diffusion coefficient and the advective wind. 

In this Letter we present for the first time a consistent 
picture based on that scenario that aims to overcome all 
of the aforementioned problems. 

The most significant achievements we present are: 

• a natural explanation to the long-standing Milagro 
anomaly; 


• a new prediction of the Galactic neutrino diffuse 
emission that is significantly larger than the one 
computed with conventional models; and 

• a possible interpretation for the hints of an excess of 
IceCube events along the Galactic plane and of the 
different neutrino slope in the northern and south¬ 
ern hemispheres. 

2. THE MODEL. 

Following dGaggero et al.l 12014! ). the starting point is 
a conventional propagation setup characterized by 5 = 
0.50. compatible with a Kraichnan spectrum of the inter¬ 
stellar turbulence within the quasi-linear theory frame¬ 
work. We will refer to thi s setup as the “KRA model” 
(see also (lEvoli et al.ll201lD l. 

The new model presented in that paper features <5 in¬ 
creasing with the galactocentric radius R (implying spa¬ 
tiallysuggested, e.g., 
in (lErlvkin fc Wolfendald 1201311 1. and hence predicts a 
hardening of CR propagated spectrum and 7 -ray emis- 
sivity in the inner Galaxy. The following explains the 
model in more detail: 

• S has the galactocentric radial dependence S(R) = 
AR + B for R < 11 kpc where A = 0.035 kpc ^ 1 
and B = 0.21 so that S(Rq) = 0.5. This behavior 
may have different physical interpretations, e.g. a 
smooth transition between a dominant parallel es¬ 
cape along the poloidal component of the regular 
Galactic magnetic field (in the inner Galaxy, where 
<5 is lower) and a perpendicular escape with respect 
to the regular field lying in the plane (in the outer 
Galaxy, where the scaling is steeper). 

• An advective wind for R < 6.5 kpc with velocity 
Vc{z)z (z is the distance from the GP) vanishing at 
z = 0 and growing as dV c /dz = 100 km s -1 kpc ^ 1 
is also included. This ingredien t is motivated 
by the X-ray ROSAT observations iSnowden et ahl 

dS) 

• The vertical dependence of the diffusion coefficient 
is taken as D(z) oc exp (z/zt)', 

• The halo size is Zt = 4 kpc for all values of R 
(this is a conventional choice widely used in the 
literature, and we checked that our results do not 
change significantly considering larger values of z t ). 

The observed 7 -ray spectra at both low and mid Galac¬ 
tic latitudes, including the Galactic center, are repro¬ 
duced by this model without spoiling local CR observ¬ 
ables: proton, antiproton and helium spectra, B/C and 
10 Be/ 9 Be ratios. Moreover, this scenario naturally ac¬ 
counts for the radial dependenc e in the CR spec trum 
found by the Fermi collaboration (lCasandiianll2015l) . We 
will refer to this model as “KRA 7 ” since it is tuned on 
gamma-ray data. 

We implement the setup with DRAGON, a numeri¬ 
cal code desi gned to compute the propagation of a ll 
CR species ( Evoli et~ahl 120081 iDRAGON-webl 1201511 . 

1 <5 is defined by D(p) oc ( p/po) S 
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While th e current version of the code shares with 
GALPROP (1GALPROP -we B I20l5h the same spallation 
cross-section routines and the gas distribution, its inno¬ 
vative structure allows us to compute CR transport in 
the general framework of position-dependent diffusion. 

Concerning the p a nd He spectral ha rdening in¬ 
ferred from PAMELA (lAdriani et al.l 12 01 111 - recently 
confirmed by AMS-02 (lAguilar et al.l I20l7j|) - and 
CREAM (I Alin et al.l I2010D data above ~ 250 GeV/n, 
we consider two alternatives. 1) Local hardening 
could ori ginate from nearby supe rnova remnants (see 
e.g. (iThoudam fc Horandelf I2013j) ): since this is a 
stochastic effect and averages out on large scales it 
amounts to not introducing any feature in the Galactic 
CR population used in this work. 2) Global hardening 
could originate from a spectral feature in the rigidity de¬ 
pendence of CR source spectra or the diffusion coefficient 
(here we only consider the former case, as both scenar¬ 
ios have the same effect on the 7 -ray diffuse emission). 
In both cases we assume that above 250 GeV/n the CR 
source spectra extend steadily up to an exponential cut¬ 
off at the energy E cut /nucleon. 

We consider two representative values of this quantity, 
namely A? cut = 5 and 50 PeV which - for the KRA 7 setup 
- match CR EAM p and He data and roughly bracket 
KAS CADE (lAntoni et al.ll2005l) and KASCADE-Grande 
data (lApell 120131) . While KASCADE proton data fa¬ 
vor the lowest cutoff, the highest is favored by the 
KASCADE-grande all-particle spectrum. A more de¬ 
tailed fit of the CR spectra in the PeV region is not jus¬ 
tified here due to the large experimental uncertainties on 
the elemental spectral shapes and normalizations. The 
consequent uncertainty on the neutrino flux should, how¬ 
ever, be well captured by our choice to consider a range 
of cutoffs. 

3. THE 7 -RAY SPECTRUM. 

As shown in dGaggero et al.l[20l4l) . the KRA 7 setup 
both in its local (KRA 7 with no hardening) and global 
realizations - provides a good fit of the 7 -ray diffuse 
emission measured by Fermi -LAT all over the sky, par¬ 
ticularly toward the inner GP region. Moreover, it ac¬ 
counts for the galactocentric radial dependence of the 
CR s pectral index fou nd by the Eerrm-LAT collabora¬ 
tion (lCasandiianll2015l) . 

Here we extend the computation performed in 
(jCaggero et al. 120141) above the TeV . 

Similar to ( Gaggero et al.l 12014 ?). we compute the 
hadronic emission integrating the expression of the 7 - 
ray emissivity along the line of sight using GammaSky, 
a ded icated code used in (lEvoli et al.ll2012l : lCirelli et al l 
l2014f ) to simulate diffuse 7 -ray maps. This package fea¬ 
tures, among other options, the gas maps included in 
the public version of GALPROP ( iGALP ROP-webl 120151 : 
iMoskalenko et ahl l2002t lAckerma nn et a l.l 120121)! We 
adopt the emissivities given in dKamae et al .112001)1) . ac¬ 
counting for the energy dependence of the pp inelastic 
cross section (significant above the TeV). We disregard 
7 -ray opacity due to the interstellar r adiation field, s ince 
it is n egligible up to a few tens of TeV (|Ahlers fc Murasel 

l20H . 

Our results are shown in Fig. [T] As mentioned in 
the above, a representative conventional model (KRA) 
cannot account for the flux measured by Milagro from 


the inner GP at 15 TeV even if accounting for the CR 
spectral hardening required to match the PAMELA and 
CREAM data. The KRA 7 setup, instead, is more suc¬ 
cessful, especially if a global hardening is assumed. This 
is a remarkable result since: 1) it supports the KRA 7 
model in a higher-energy regime; 2) it provides the first 
consistent interpretation of Milagro and Fermi -LAT re¬ 
sults (an optimized mod el was p r opos ed to account for 
the EGRET GeV excess [Strong et all (120041) . and came 
out to reproduce Milagro results as well, but was su b- 
sequently excluded by Ferrm-LAT lAbdo et~aT1 (j2009D ). 
and 3) it reinforces the arguments in favor of a non-local 
origin of the hardening in the CR spectra above 250 GeV. 

Interestingly, the KRA 7 model also reproduces the 
high-energy diffuse 7 -ray spectrum measured by H.E.S.S. 
in the Galactic ridge region (|Z| < 0.8°, | 6 | < 0.3°) in 
terms of CR scattering with the dense gas in the central 
molecular zone without the nee d to invoke the contri- 
bution of sources to that region (lAharonian et all 120061) 
and without further tuning (see dGaggero et al.ll2015l) for 
more details). Although this is a very small region with 
respect to the regions considered in this paper, this re¬ 
sult may be interpreted as a valuable check of our model 
in a region not covered by Milagro 

Moreover, our KRA 7 model is also compatible with 
ARGO-YBJ results in the window 65° < l < 85° and 
| 6 | < 5°; both the KRA and the KRA 7 are consistent 
with CASA-MIA measure ments at high Galactic longi¬ 
tudes dBorione et alJflflflM ). 

4. THE NEUTRINO EMISSION. 

The results discussed above clearly show that the 
hadronic emission computed with the KRA 7 setup above 
the TeV is significantly stronger than the conventional 
model predictions. In this section we show the relevant 
consequences concerning the Galactic neutrino emission. 

We first compute the v e and production spec¬ 
tra: for both f la vors we use the emissivities provided 
in (|Kamae et al.l I2006D (well tuned on accelerator and 

CR data) for projectile energies below ~_500 TeV, 

while we adopt the emissivities provided in ( iKelncr et al.l 
i 200(1 ) that are above that energy range. Then we ac¬ 
count for neutrino oscillations: their effect is to almost 
equally redistribute the com position among the three fla¬ 
vors dCavasinni et al.l [20061) . We only consider proton 
and helium CRs/gas - as for 7 -rays - since heavier nu¬ 
clear species give a negligibl e contribution in t he energy 
range we cover in this work dKachelriess fc Ostapchenkol 
12014) . 

Because neutrinos in the Galactic emission are ex¬ 
pected to be maximal in the inner Galactic plane re¬ 
gion, we first present our results for the window |Z| < 30° 
and | 6 | < 4°. For this reg ion the ANTARES collabo¬ 
ration (| Aguilar et al.l 1201 11 recently released an upper 
limit on the muon neutrino flux based on the result 
of an unblinding analysis regarding the events collected 
betw een 2007 and 2013 in the energy range [3 4- 300] 
TeV (lEnseo fc ANTARESII 20 T 5 I) . 

In Fig. [2] we compare the flux computed with the 
KRA and KRA 7 setups with the flux of the experimen¬ 
tal constraint. First of all we notice the large enhance¬ 
ment (almost a factor of 5 at 100 TeV) obtained with 
the KRA 7 model with respect to the conventional sce¬ 
nario. Indeed, while - in agreement with previous re- 
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Figure 2. Solid and dashed red (blue) lines: expected neu¬ 
trino spectra (all flavors, both neutrinos and antineutrinos) in the 
inner Galactic plane region computed for the conventional KRA 
(the novel KRA 7 ) models for two different cutoff values. We also 
show the maximal flux, estimated c onsidering th ree years of Ice- 
Cube HESE events as described in (Spurio 2014 ), the con straint 
from the ANTARES experiment (IFusco fe ANTARESII20151) (1500 
days of experiment livetime between 2007 and 2013) as well as 
the deduce d sensitivity of the future Mediterranean observatory 
KM3NeT JPiattelli &s KM3NeT1 120151) with four years (~ 1500 
days) of livetime. 


Figure 3. Full-sky neutrino spectrum (all flavors, both neutri¬ 
nos and antineutrinos) predicted by the KRA 7 and KRA mod¬ 
els (with global CR hardening), adopting two different choices for 
the CR high-energy cutoff. We also plot the combination of the 
Galactic (KRA 7 ) and a benchmark EG spectrum. The EG flux 
is consistent with that inf erred from the IceCub e collaboration in 
the northern hemisphere (lAartsen et al.lf2015bl) . The models are 
compared with the 68 % confidence region for the IceCube astro- 
physical neutrino flux obtained with a maximum-likelihood (yellow 
region) JAartsen et al.ll2015al) and the three years HESE (green 
points) <|Aartsen et al.1120141) . 


suits - we find that the flux corresponding to the KRA 
model may require lon g times of observation even b y the 
KM3NeT observatory (jAdrian-Martmez et al.ll2013[ h our 
prediction for the KRA 7 model is instead well above the 
sensitivity reachable by that experiment in four years and 
it is almost within the ANTARES observation capabili¬ 
ties. 

Interestingly, our result is in good agreement with the 
maximal flux inferred from the fraction of IceCube HESE 
events compatible with that region (see Fig. We no¬ 
tice that in that region the expected EG contribution, as 
constrained from the muon neutrino flux in the northern 
hemisphere (see below) gives a subdominant contribution 
with respect to that computed with the KRA 7 model. 
Therefore the possible detection of a signal in that sky 
window would be a smoking gun for the presence of such 
Galactic emission. 

IceCube should also have the potential to detect that 
emission on a larger region. In t his context, we also 
note t hat an independent analysis dNeronov fc Semikozl 
I2015al) already found a significant hint of an excess in 
the 4-year HESE sample (lAartsen et al.ll2015dll along the 
Galactic plane. 

We now turn our attention to the recently published 
IceCube results, both concerning the full-sky and the 
northern and southern hemispheres separately. 

In Fig. [3] we represent the full-sky total neutrino spec¬ 
trum (all flavors, including antiparticles) computed for 
the KRA 7 and KRA models, with global CR hardening, 
and compare it to the IceCube results. 

Our prediction for the convent ional setup (KR A 
model) is in good agreement with (jAhlers ct all 120151) : 


In that work, the benchmark Galactic model accounts 
for 8% of the flux measured by IceCube above 60 TeV, 
for a CR spectrum similar to the one used here above 50 
PeV. 

On the other hand, the KRA 7 predicts a ~ 2 times 
larger full-sky flux above 10 TeV: the model prediction 
is therefore only ~ 4 times smaller than the best fit of 
the astrophysical flux measured by IceCube on the whole 
sky. 

We r emark that another analysis dNeronov fe Semikozl 
I2015IJ) . based on an extrapolation of Fermi -LAT data, 
points toward a non-negligible Galactic contribution to 
the full-sky neutrino flux due to a hard diffuse CR spec¬ 
trum. In that scenario the (softer) locally observed CR 
spectrum may get a major contribution from one or more 
local sources: this interpretation still has to be validated 
against Fermi -LAT data, while our model is based on 
those measurements. 

Setting a threshold energy at 60 TeV and convolving 
the KRA 7 spectrum (with E „ ltl = 50 PeV) with th e Ice- 
Cube HESE effective areas (lAartsen et al.|[2013alh the 
expected number of neutrino events in three years of 
IceCube observation represe nts ~ 15% of the published 
sample dAartsen et al.|[2bl4f) . These rates are well above 
those expected due to atmospheric muons and atmo¬ 
spheric neutrinos and confirm the spectral comparison 
between KRA 7 and IceCube data. 

Clearly, another component - most likely of extragalac- 
tic (EG) origin - needs to be invoked in order to account 
for all of the IceCube events. 

Here we assume this EG component to be isotropic 
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and use the astrophysical muon neutrino IceCube mea¬ 
su rem ents from the northern hemisphere (lAartsen et al.l 
l2015bT) - where the Galactic emission is only ~ 1/10 
of the total flux - to probe its spectral properties. Al¬ 
though the northern spectral slope is statistically com¬ 
patible with the full-sky one, given the hint of a steeper 
spectrum in the southern hemisphere, it is interesting to 
check if the combination of our Galactic prediction and 
the EG flux inferred from the aforementioned muon neu¬ 
trino measurement provide a better agreement with the 
data. 

For illustrative purposes, in Fig. [3] we show the effect 
of adding an isotropic EG emission to the Galactic neu¬ 
trino emission computed with the KRA 7 model, with a 
spectrum given by the IceCube best fit of c E > /?° rth , mul¬ 
tiplied by three to account for all flavors. The nature 
of such emission is still under deba te: as pointed out i n 
(iGliisenkamp fc IceCubef2015 ) and (|Bechtol et al.ll20T5l b 
neither blazars nor star-forming galaxies can provide 
more than a subdominant contribution, given the con¬ 
straints imposed by the gamma-ray extragalactic back¬ 
ground inferred from Fermi -LAT data. The plot clearly 
shows how the KRA 7 helps to improve the fit in the low- 
energy part of the IceCube spectrum. 

We also checked that the neutrino flux computed with 
the KRA 7 model for | 6 | < 7.5° is in rather good agree¬ 
ment with that inferred from IceCube HESE analysis if 
the EG emission, as estimated above, is accounted for. 
A dedicated analysis will be performed in a forthcoming 
work. 

5. CONCLUSIONS. 

In this Letter we connected 7 -ray GeV and TeV mea¬ 
surements in a unified scenario, together with the re¬ 
cently released IceCube neutrino data, providing a con¬ 
sis tent picture based on a CR transport model proposed 
in flGaggero et al.ll2014f) . The model features a variation 
of the diffusion coefficient rigidity scaling S with galacto- 
centric radius. The variation was suggested by a spectral 
anomaly found in the Fermi -LAT 7 -ray data, and turned 
out to be compatible with both 7 -ray spectra at low and 
intermediate Galactic latitude and local CR observables. 

In this work we showed that our picture sheds new 
light on high-energy gamma-ray and neutrino recent ob¬ 
servations. In particular, it provides a novel natural ex¬ 
planation for the anomalous 7 -ray flux measured by the 
Milagro observatory from the inner GP region at 15 TeV; 
moreover, it appears to be compatible with the H.E.S.S. 
spectrum in the Galactic ridge region. 

Remarkably, our model also provides a different inter¬ 
pretation of the full-sky neutrino spectrum measured by 
IceCube with respect to the standard lore, since it pre¬ 
dicts a larger contribution of the Galactic neutrinos to 
the total flux, compared to conventional models. 

These predictions will be testable in the near future 
by neutrino observatories such as ANTARES, KM3NeT, 
and IceCube itself via dedicated analyses that are fo¬ 
cused on the Galactic plane, and also by analyzing the 
different spectral slopes in the northern and southern 
hemispheres. A hint of a softer slope in the northern 
hemisphere is already present, and appears to be com¬ 
patible with our picture. 

A physical interpretation of our model most likely re¬ 
quires either abandoning the isotropic diffusion scenario 


generally adopted to treat CR propagation, or consider¬ 
ing different turbulence regimes in different regions of the 
Galaxy: a quantitative modeling of those phenomena is 
far beyond the scope of our phenomenological work. 
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